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In the last 20 years, heavy-ion collisions have been a unique way to study the 
hadronic matter in the laboratory. Its phase diagram remains unknown, although 
many experimental and theoretical studies have been undertaken in the last decades. 
After the initial experiences accelerating heavy nuclei onto fixed targets at the AGS 
_~ (BNL, USA) and the SPS (CERN, Switzerland), the Relativistic Heavy Ion Collider 

(RHIC) at BNL was the first ever built heavy-ion collider. RHIC delivered its first 
collisions in June 2000 boosting the heavy-ion community. Impressive amount of 
experimental results has been obtained by the four major experiments at RHIC: 
PHENIX, STAR, PHOBOS and BRAHMS. In November 2010, the Large Hadron 
Collider (LHC) at CERN delivered lead-lead collisions at unprecedented center-of- 
mass energies, 14 times larger than that at RHIC. The three major experiments, 
ALICE, ATLAS and CMS, have already obtained many intriguing results. Needless 
to say that the heavy-ion programs at RHIC and LHC promise fascinating and 
?j exciting results in the next decade. 

The first part of the lectures will be devoted to introduce briefly the QCD descrip- 
tion of the strong interaction (as part of the Standard Model of Particle Physics) 
and to remind some basic concepts on phase transitions and on the phase diagram 
of matter. 

In the second part, I will focus on the properties of matter at energy densities above 
~1 GeV/fm 3 . A historical approach will be adopted, starting with the notion of 
(NJ limiting temperature of matter introduced by Hagedorn in the 60's and the discovery 

of the QCD asymptotic freedom in the 70's. The role played by the chiral symmetry 
breaking and restoration in the QCD phase transition will be discussed, supported by 
an analogy with the ferromagnetic transition. The phase diagram of hadronic matter, 
conceived as nowadays, will be shown together with the most important predictions 
of lattice QCD calculations at finite temperature. Finally, the properties of an 
academic non-interacting ultra-relativistic QGP and its thermal radiation will be 
t> deduced. The dissociation of the heavy quarkonium due to the color-screening of the 

heavy-quark potential will be described, based on a QED analogy. The energy-loss 
phenomenology of ideal long-living partons traversing the QGP, will be reminded. 

In the third part, the heavy-ion collisions at ultra-relativistic energies will be 
proposed as a unique experimental method to study QGP in the laboratory, as 
suggested by the Bjorken model. The main experimental facilities in the world will 
be described, namely the CERN and BNL accelerator complexes. The main probes 
for characterizing the QGP in heavy-ion experiments, followed by a brief description 
of the main heavy-ion experiments located at these facilities will be shown. 

In the last part of these lectures, I will present my biased review of the numerous 
experimental results obtained in the last decade at RHIC which lead to the concept 
of strong interacting QGP, and the first results obtained at LHC with the 2010 and 
2011 PbPb runs. Finally, the last section is devoted to refer to other lectures about 
quark gluon plasma and heavy ion physics. 
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I. INTRODUCTION 

The strong interaction, described by quantum chromodynamics (QCD), is the dominant 
interaction in the subatomic world. The main properties of the strong interaction are: 

• the strength constant a s at low energies is large, and as a consequence, quantitative 
calculations based on a perturbative sum of Feynman diagrams, faiQ 

• gluons, g (the intermediate boson of strong interaction) are coloured (colour is the 
charge of the strong interaction). For this reason, the QCD becomes a complex quan- 
tum field theory (QFT), which belongs to the class of non-abelian QFT. 

The quarks u, d and s, also called light quarks, exhibit small masses and therefore the 
most important parameter of QCD is indeed a s . However a s can only be determined in 
the high energy domain, since its experimental determination at low energy is difficult due 
to non-perturbative effects. One of the major experimental observations that QCD should 
explain, is the confinement of quarks and gluons. Coloured free particles do not exist, and 
thus quarks and gluons seem to be confined inside colourless particles called hadrons. The 
confinement property is not fully understood, despite the fact that the quark model describes 
qualitatively the hadron properties (mesons are bound states of a quark and antiquark and 
baryons are bound states of 3 quarks). Today, the best ab-initio quantitative calculations 
can be performed via lattice calculations of QCD. One should note that the origin of the 
hadron mass is the strong interaction, since light quark masses only represent less than 10% 
of the total hadron mass. As Frank Wilczek (Nobel prize in 2004) expressed in Physics Today 
in November 1999 "According to quantum chromo- dynamics field theory, it is precisely its 
color field energy that mostly make us weigth. It thus provides, quite literally, mass without 
mass". In this respect the Higgs boson (strictly speaking the Brout-Englert-Higgs boson or 
BEH boson) only explains about 1% of the total mass of the proton and neutron which are 
the main massive constituents of ordinary matter. 

A. Asymptotic freedom 

The vacuum polarization of QCD |Politzer 73[ IGross 73] exhibits a singular behaviour 
due to the anti-screening effect of virtual gluon pair production (remember that gluons are 
colored bosons and the gluon vertex does exist in QCD). Indeed the gluon anti-screening is 
stro nger that the screening effect of virtual quark pair production (see Fig. [I]). In QCD one 
gets [Griffiths 87] : 

«(k 2 |) = (2) asW (1) 

l + ^(lln-2/)ln(|gW) 

where n represent the number of colors and / the number of quark flavors. In nature, 
lira > 2/, in consequence, the strength of the strong interaction a s decreases at small 
distances (or high energies) . This phenomenon is called the asymptotic freedom of QCD. 
The discovery of the asymptotic freedom was awarded with the Nobel Prize in 200^ 
At the scale of the Z boson mass (|g 2 | ~ Mz), oc s has been measured via many different 



physics channels and the current world average is 0.1184±0.0007 Beringer 12 . Perturbative 
QCD calculations are then fully valid to describe the strong interaction at high energies. 
This is one of the major successes of QCD theory. 



1 Roughly, the strength constant a s can be estimated from the hadron bound state properties. Using the 
Bohr radius expression of the hydrogen atom (please forgive me this unacceptable assumption) rs — 
2 /(ma), one can estimate a s ~ 20, considering tb = 1 fin and m q = 1 MeV. 



http : //www . nobelprize . org/nobel_prizes/physics/laureates/2004/ 





FIG. 1. Feynman diagrams at first order, of the vacuum polarization in QCD: a) screening and 
b) anti- screening. In the case of QED, anti- screening does not exist since photons are not charged 
particles. 

In equation (fTl), the parameter /i 2 is imperative since perturbative QCD calculations 
cannot be performed in the domain where a s ^ 1 (when \q 2 \ < 1 GeV 2 ). The expression (fTl) 
can be rewritten as a function of a parameter Aqcd'- 



«(k 2 l) 
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Aqcd is 213 ± 8 MeV for 5 flavours |Beringer 12| . However, the quantity Aqcd is not well 



defined. Therefore it has become standard practice to quote the value of a s at a given scale 
(typically Mz) rather than to quote a value for Aqcd |Beringer 12 . 



B. Lattice QCD calculations 



As we have previously seen, the asymptotic freedom is the main success of QCD and it 
allows for the experimental test of QCD via the study of high energy processes. I would say 
that the description of i) the evolution of the parton distribution functions at low Bjorken x 
values, ii) the production of jets in elementary collisions and iii) the properties of bottomo- 
nium bound states , represent beautiful examples of the success of QCD predictions in its 
perturbative domain. 

Nevertheless QCD should a priori explain many other phenomena where perturbative 
calculations are forbidden. Fundamental questions such as the coupling constant value, 
QCD vacuum structure, hadron masses, hadron structure, nuclear properties etc ... should 
be explained by the QCD theory. Even today, the test of QCD has not been possible in 
many of these domains. In order to palliate this, many effective models have been developed 
in the domain of hadronic and nuclear physics. 

Lattice calculation for gauge theories is the most promising non-perturbative technique 
to solve QCD equations. The space-time continuum is discretized in a finite number of 
points where the equations of the theory can be solved. In the last decade, the impres- 
sive development of computing hardware and the optimization of software algo rithms have 
allowed lattice QCD calculations to become a competitive tool [Wilczek "03] . Nowadays 
computing facilities for lattice QCD calculations are a crucial component of this research, 
at the same level as accelerators, detectors or computing centers for data analysis and data 
storage. Today lattice QCD provides the most precise computation of the a s constant, is 
able to extract the mass of the quarks and to predict the mass of most of the hadrons, and 
it makes excellent predictions of the exotic structure of new bound states of heavy quarks. 
In particular, we will see later in more detail, how lattice QCD at finite temperature allows 
to study the hadronic matter phase diagram. 



C. A description of the hadronic matter phase diagram 

By hadronic matter I mean that in which the strong interaction is the main interaction 
between elementary constituents, that provide the proper degrees of freedom of the matter. 
At temperatures above 10 9 K (1 MeV) and/or pressures above 10 32 Pa (1 MeV/fm 3 ), the 
strong interaction is expected to be the dominant interaction between the constituents of 
matter. At low temperatures and a pressure above 1 MeV/fm 3 the matter can be described 
as a degenerated gas of neutronq^j Such a state, which is very close to the atomic nucleus 
structure, should exist in the neutron stars. In these stellar objects, a mass slightly larger 
that the sun mass, is confined in a ten kilometer radius sphere, and densities as high as 10 17 
kg/m 3 are reached. For higher pressures, above 10 35 Pa (1 GeV/fm 3 ), the repulsive force of 
the degenerate gas of neutrons cannot compensate the pressure, and matter is expected to 
become a low temperature gas of quarks which are not any more confined inside hadrons. 
In this exotic state of matt er, quark-quar k Cooper pairs might exist creating a kind of 
color superconductor matter |Rajagopal 00| . On the other hand, the neutron matter should 
become a gas of nucleons, if it is heated to temperatures of several MeV. Indeed the nucleon- 
nucleon potential has some similitudes with the Van der Waals force between molecules. For 
this reason, it is expected that the neutron matter evaporates into a gas of nucleons at a 
temperature of about 10 MeVQ like the liquid-gas phase transition in ordinary matter. 

At very high temperatures and pressures, the nucleon gas (that has become a hadron 
gas at temperature above 100 MeV) could go through a transition to a deconfined state of 
matter. This is expected due to the vacuum polarization at the origin of the asymptotic 
freedom of QCD. Therefore the strength of the strong force decreases at high temperature. 
The deconfined state of matter, in analogy with the electromagnetic plasma where ions and 
electrons are dissociated, has been called Quark Gluon Plasma (QGPFJ The transition to 
QGP takes place at temperatures about 200 MeV(~2-10 12 K), when quarks and gluons are 
not confined in colorless particles and they become the pertinent degrees of freedom of the 
system. Other properties of QCD also predict the reason that a phase transition should occur 
at high temperature. In quantum field theories, the symmetries of the Lagrangian can be 
spontaneously broken at low energies or temperatures. In the case of QCD, the spontaneous 
breaking of the chiral symmetry takes place at low temperature. The restoration of the 
chiral symmetry at high tem peratures becomes a sufficient condition for the existence of a 
phase transition Smilga 03| . 



Finally, for temperatures above 10 16 K, it is hard to know what would be the structure of 
matter. Some authors have speculated about new phenomena like formation of microscopic 
black holes or unification of interactions, etc ... that could appear. Exotic i deas like th e 
formation of a superstring gas have been proposed for temperatures of 10 32 K |Bowick 85] . 

In Fig. |2j the lay-out of the phase diagram of matter is presented. We clearly distinguish 
two regions, one for temperatures below 10 9 K and pressures below 10 30 Pa, where the 
electromagnetic interaction between atoms (or ions) provides the degrees of freedom of 
matter, and a second region, for temperatures above 10 9 K and/or pressures above 10 32 Pa, 
where the strong interaction between nucleons, hadrons or quarks dominates and provides 
the right degrees of freedom. 



3 Actually, it is a degenerate gas of baryons and electrons which is a more stable system that a pure neutron 
gas. 

4 Naively, this value can be accepted since the bound energy of nucleus in the nuclei saturates to a value of 8 
MeV. The liquid-gas transition has been studied in heavy ion collisions at intermediate energies, although 
it is not discussed in these lectures. 

5 Professor E. Shuryak proposed the name Quark Gluon Plasma in the 80's Shuryak 78 . 
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FIG. 2. Phase diagram of matter in the pressure versus temperature plane for a non zero baryonic 
potential [Martinez 06]. 



II. THE QUARK GLUON PLASMA 



A. Limiting temperature of matter 

Curiously, the first prediction of a critical behaviour of hadronic matter at high tempera- 
ture was obta ined before the formulation of QCD and the discovery of partons |Hagedorn 65 



Hagedorn 84| . At mid 60's, Hagedorn was interested in properties of the hadron gas. He 



predicted in a phenomenological and or iginal manner , that should exist a critical behaviour 
of the hadron gas at high temperature jHagedorn 65] . Hagedorn interpreted this criticality 
as the existence of a maximum temperature of matter, that was called Hagedorn's temper- 
ature Th- In order to study the hadron gas, one has to consider all the zoology of hadron 
particles. Today, more than 2000 hadron species have been discovered (see Fig. [3] left). 
Hagedorn studied the number of hadron species as a function of their mass. He observed an 
exponential dependence and the following function was used to describe the experimental 
data: 



p(m) 



A 



m 2 + [500MeV J 



exp (m/T H ) 



(3) 



where p(m) is the density of hadron species per mass unit and T# is a parameter. From the 
experimental data, one obtains that the parameter Th is close to the mas s of the pion, ~ 180 
MeV, when all the known baryon and meson resonances are considered [Broniowski 04] . It 
turns out that such a dependence of the density p{m) will induce divergences of the partition 
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FIG. 3. Left: number of hadron species as a function of their mass \Broniowski 0$ . Right: first 
phase diagram of hadronic matter fCabibbo 75f . ps is the baryonic density, T the temperature, (I) 
is the confined phase and (II) is the deconfined phase. 



function that describes the statistical properties of a hadron gas, if the temperature of matter 
reaches values above the Hagedorn parameter T#. In consequence, Th was interpreted as a 
limiting temperature of matter. Somehow, any additional energy supplied to the system at 
the Hagerdorn temperature, would be used to create new hadron species. 

We know that hadrons are not point-like particles, and their typical size is around 1 fm 
sphere radius. Indeed, when one gets closer to the T H temperature, the hadron density 
increases (remember that the energy density of an ideal ultra-relativistic gas increase as T 4 ) 
and values about 1 hadron per fm 3 are reached. Under these conditions, hadrons overlap 
with each other and considering hadrons as point-like particles (which means that their size 
is small with respect to its mean free path) becomes a wrong hypothesis and invalidate Hage- 
dorn conclusions. Therefore one has to understand first the internal structure of hadrons, 
since it is going to provide the new degree of freedom of the system when T> Th- Only 
QCD was able to answer this question several years later. 



B. Deconfined state of matter 



After the discovery of the asymptotic freedom |Politzer 731 1 Gross 73] , the existence of 
a deconfi ned state of quarks and gluons was predicted at high temperature and/or high 
pressures |Collins 75} ICabibbo 75j . A fir st pioneer ph ase diagram of hadronic matter was 
imagined (see Fig. [3] right from reference [Cabibbo 75j ). At sufficiently high temperatures, 
quarks and gluons interact weakly and the system will behave as an ideal ultra-relativistic 
gas. The degrees of freedom will be then determined by the flavor numbers, spin states, 
color and charge st ates of the q uarks and gluons. The deconfined state was called later 
quark gluon plasma |Shuryak~78 . The word plasma is used to describe the state of matter 
when ions and electrons are dissociated in atoms. There is then an analogy when colourless 
particle dissociate to create deconfined matter. One open question after the discovery of 
the asymptotic freedom, concerned the properties of the transition from the hadron gas 
to the QGP: does it take place smoothly or via a phase transition and exhibiting critical 
behaviours? As a matter of fact, the transition from gas to elec tronic plasma takes place 
smoothly in the temperature range 10000 to 50000 K [Stocker T)9] and no critical behaviour 
is therefore observed. The question whether the QGP phase transition exists, is, of course, 
a very deep question and the intrinsic symmetries of the QCD could give us the answer. 
Indeed the chiral symmetry of the massless quark QCD Lagrangian is spontaneously broken 



at low temperature and this symmetry should be restored at high temperatures. A symmetry 
restoration represents a valid condition to predict the existence of a QCD phase transition. It 
remained however an open question if the chiral symmetry transition and the deconfmement 
transition are or not the same one. Only lattice calculations have been able to provide an 
answer to this question as we will see later. 

C. The spontaneous break-up of chiral symmetry in QCD 
A simplified Lagrangian of 3 quark flavours / (u, d ,s) can be written as [Schaefer 05] : 



N f 



4 
/ 



(4) 



where Nf = 3 and the coupling gluon field tensor is defined as: 



G% = d,Al - d v Al + gr hc A\A^ (5) 

and the covariant derivate of the quark field as: 

m = ^Ud^ + gA a ^\il>. (6) 

Under these conditions, the previous Lagrangian exhibits a flavour symmetry since the 
quark interaction does not depend on the quark flavour. This is indeed always the case 
if the masses of the quarks are identical. The direct consequence of this is the symmetry 
under isospin transformations, that it is observed in the hadron properties. In addition, for 
massless quarks, the QCD Lagrangian exhibits the chiral symmetry^) The quark fields can 
be decomposed in left-hand and right-hand quarks fields [Halzen 84] : 

Kr=^(i± 75 )v>. ( 7 ) 

As a consequence, the QCD Lagrangian is invariant under helicity and flavour transforma- 
tions. This symmetry is represented as the SU(3)l x SU(3)r symmetry of QCD. One of 
the consequences of this symmetry is that the associated parameter, called condensate (qq) 
should be zero. 

Nevertheless, the condensa te (qq) is n ot zero and the existence of the pion is a clear 
confirmation of this statement |Knecht 98] . This is what it is called the spontaneous breaking 
of the SU(3)l x SU(3)r chiral symmetry of QCD. The word spontaneous reminds us that the 
symmetry is respected by the QCD Lagrangian but broken by their states at low energies. 
At high energies the symmetry should be restored. 

The spontaneous breaking of a symmetry is a phenomenon that is allowed in quantum 
field theories, where the structure of the vacuum plays a major role. In quantum mechanics, 
the eigenstates that respect the symmetry of the Hamiltonian, can always be found. In 
classical mechanics the following analogy of the spontaneous symmetry breaking can be 
found. Let's assume a ring in the earth gravitational field, that rotates along its vertical 
symmetry axis with an angular speed u. There is a small solid ball with a hole in a manner 
that can move freely along the ring (see Fig. 41). In this example, the system exhibits a 
left-right symmetry which is spontaneously broken by the small ball at low internal energy. 



Chiral from hand in Greek. 




FIG. 4. Classical analogy of spontaneous symmetry breaking. The ball is holed and can move 
freely along the ring, which rotates with an angular speed w, in the earth gravitational field with 
respect to its vertical symmetry axis. 



Indeed, due to the centrifugal force, the small ball has to choose the left or the right side of 
the ring as its equilibrium position. If some internal energy is given to the ball, it will start 
to oscillate around its equilibrium position. The amplitude of the oscillation will increase 
with the internal energy of the ball. Above a certain energy threshold, the ball will have 
enough internal energy to reach the other side of the ring and it will then move in both 
sides. When this occurs, one can say that the left-right symmetry of the system has been 
restored. 

The spo ntaneous breaking of the chiral symmetry is one of the predictions of QCD 
|Knecht 98] . and, in this way, QCD is able to predict the existence of the Goldstone bosons: 
the pions, kaons and eta mesons and to explain their small interaction cross-sections. As 
in the classical analogy, the chiral symmetry of the QCD is restored at high energies (or 
high temperatures) and remember that a restoration of the symmetry represents a sufficient 
condition for the existence of a QCD phase tran sition. A n analogy with the ferromagnetic 
phase transition can be made (see table III) |Schaefer 05] . In fact, the ferromagnetic phase 
transition can be associated to the spontaneous breaking of the isotropy symmetry. At high 
energy the ferromagnetic system is invariant under rotation transformation, since there is 
not any privileged direction of the space. Nevertheless at low temperatures, the thermal ag- 
itation cannot avoid that the microscopic magnetic moments of the elementary constituents 
align, causing a macroscopic magnetisation of the system. Therefore the isotropy symmetry 
is spontaneously broken at low temperatures, and this is a sufficient condition to predict 
that there is a phase transition during the generation of the macroscopic magnetisation of 

the system. In the ferromagnetic case, the magnetisation M is the order parameter of the 
transition, which is the equivalent of the quark condensate (qq) in the chiral transition in 

QCD. The non-zero M, allows for the existence of spin waves, and the Goldstone bosons 
(pions, kaons and eta's) are their analogous. Finally, isotropy symmetry can be explicitly 
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Transition 


Chiral 


Ferromagnetic 


Spontaneous breaking 
Order parameter 

States 
Explicit breaking 


SU(3) L x SU(3) R 

Condensate < qq > 

Goldstone bosons n, K, ... 

Quark masses m q ^ 


Isotropy 0(4) 

Magnetisation M 

Spin waves 

External magnetic field 



TABLE I. Analogy between the chiral and the ferromagnetic phase transition [Schaefer 05 1 . 

broken via an external magnetic field. The equivalent of the non-zero external magnetic field 
would be the non-zero masses of the quarks, which explicitly breaks the chiral symmetry of 
the QCD Lagrangian. 

We have seen that a spontaneous breaking of the chiral symmetry explains why there 
should be a phase transition of the hadronic matter. We can now wonder if such a tran- 
sition is associated to the process of deconfmement of quarks and gluons leading to the 
formation of a quark gluon plasma. One could imagine that there are, indeed, two different 
phase transitions, a chiral transition and a deconfmement one that occur at different critical 
temperatures. In the next section, lattice QCD calculations will be presented since this is 
the only way to answer this question today. 

It should be noted that we have assumed a QCD Lagrangian with massless u, d and s 
quarks. This is indeed a good approximation, since the masses of the, so called light quarks, 
are small compared to Aqcd but they are not zero: m u = 2.3 ± 0.5 MeV, nid = 4.8^o' 3 
MeV and m s = 95 ± 5 MeV |Beringer 12| . In this respect the chiral symmetry is indeed 
explicitly broken by the QCD Lagrangian. Above we have assumed that if the masses are 
small compared to Aqcd this chiral symmetry should remain a good symmetry of QCD. 
However this may be a wrong assumption, in particular for the strange quark. Indeed, it is 
an open question what would be the masses value thresholds causing the damp out of the 
criticalness of the chiral transition. Above such mass thresholds, the chiral transition would 
become a cross-over and no critical behaviour would be observed in the transition. Once 
more, the lattice QCD calculations will be a unique method to study this question. 

Finally, there is a new symmetry of the QCD Lagrangian in the limit of quarks masses 
m q — > oo. The order parameter of this symmetry is called t he Polyakov line (P) which is 
directly associated to the process of deconfmement if (P)=0 [Schaefer 05] . 



D. Some results from lattice QCD calculations at finite temperature 



Today, lattice QCD calculation is a unique method to test QCD in the non-perturbative 
domain. In the last decades, many progresses have been achieved on the algorithms and 
on the computing performances. Lattice QCD allows for non-perturbative calculations with 
high reliability. 

In particular, lattice QCD should allow to study the properties of the Universe be- 
tween few ns and few /xs after the Big Bang (temperatures around 100-1 000 MeV) and 
to study hadronic matter in the core of the neutron stars |Petreczky 12| . For massless 
quarks, these calculations show a transition at baryonic potential /i# = 0, as expected from 
the spontaneous breaking of the chiral symmetry in QCD. The critical tempe rature would 
be T = 173 ± 15 MeV and the critical energy density e = 0.7 ± 0.3 GeV/fm 3 [Karsch 02aj 
(see Fig. [5]). It is also observed that above the critical temperature, the energy density is 
indeed proportional to T 4 , as for an ideal ultra-relativistic gas, but the proportionality factor 
(Stefan-Boltzmann constant) is about 20% smaller than the expected value for an ideal gas 
of gluons and massless u, d and s quarks. Perturbative calculations at higher temperatures 
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FIG. 5. Dependence of the energy density as a function of the temperature of the hadronic matter 
at null baryonic potential given by lattice QCD calculations at finite temperature. The calculations 
are performed for two massless quarks, three massless quarks and two massless quark and one (s) 
with its real mass. A transition is observed at a temperature of about 173 MeV and energy density 
of 0. 7 Ge V/frr? . For the calculations with a real s mass, the transition is faded away IKarsch 02a\[ 



are able to explain the evolution of this factor for T > 2T C [Blaizot 99] . 

The lattice QCD calculations show that for massive quarks, the phase transition could 
fade away, it would become a cross-over and no criticalness would be observed. The crit- 
icalness of the transition has been studied as a function of the quark masses (see Fig. |6|). 
In the calculations presented here, the u and d masses are considered to be identical and 
fiB = 0. It is observed that for both low and large masses, a 1st order phase transition 
is predicted. The cross-over transition occurs for intermediate quark masses. A 2nd order 
phase transition occurs in the border line between 1st order and cross-over areas. Today 
there is some consensus to believe t hat for the physical quark masses and /xs=0 there is not 
a phase transition but a cross-over [Ka rsch 02a| IKarsch 02 bp] 

The QCD lattice calculations with physical quark masses, have determined critical tem- 
peratures between 150-200 MeV. There has been some confusion about the exact critical 
temperature of the transition in the last years. The outcome was that the evaluation of the 
transition temperature, which is not a well defined parameter for a cross-over transition, 
would depend on the method used for its determination. Calculations based on chiral order 
parameter show a cross-over transition for T~155 MeV. On the other hand, the behaviour 
of the Polyakov loop suggests that colo ur screening s ets in at temperatures that are higher 
than the chiral transition temperature |Petreczky 12| . 

Finally, lattice QCD calculations have studied the order parameters of the chiral and 
deconfinement transitions (see Fig. [7]) showing that, a priori, both transitions occur at the 

7 Note that more recent references on this subject exist and they are not referenced in this lecture. 
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3 quark flavours, [ib = and assuming the mass of the u and d quarks are identical and a strange 
quark mass, m s fKarsch 02b\ \Karsch 02a\/ . 



same critical temperature. Therefore, this suggests that both transitions would be indeed 
the same transition. However, the interplay between chiral and deconfmement aspects of the 
transition appears to be more complicated than earlier lattice studies suggested. It seems 
that there is no transition temperature that can be associated with the deco nfining aspects 
of the transition for physical values of the light quark masses |Petreczky 12 



In the last decade, a lot of effort has been done to perform calculations at fiB ^ 0. 
These calculations show that there would be a critical point at //# ~ 0.75Mjv (Mat is the 
nucleon mass) where the cross-over becomes a 2nd order phase transition, and beyond it, the 
transition bec omes a 1st order phase transition between the gas of hadrons and the quark 
gluon plasma |Fodor 03] . In addition, other calculations have predicted a transition to a 
colour superconductor matter at high values of \xb (see the lay-out of the hadronic matter 
phase diagram in Fig. pj. 



E. Properties of a QGP in the ultra-relativistic limit 

At temperatures Aqcd < T < charm mass, m c and assuming that the strong interaction 
strength becomes very small, the QGP would behave as an ideal gas. Strictly speaking 
this gas will be constituted by all the elementary particles {m <T)\ leptons (electrons and 
muons), bosons (photons and gluons) and light quarks u, d and s|jand their corresponding 
antiparticles. This gas will have similar properties as the black body radiation and the 
equation of state is given by e = 3p, where e is the energy density and p is the pressure. The 
energy density will depend on the temperature following the Stefan-Boltzmann law oc T 4 . 



We assume a temperature lower than the mass of the charm quark and in this manner one could neglect 
the thermal production of heavy quarks, tau leptons and weak bosons 
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FIG. 7. Critical behaviour for massless quarks and ^b = of the order parameters of the 
deconfinement (left plot) and of the chiral (right plot) transitions as predicted by lattice QCD 
calculations. The order parameters are the Polyakov susceptibility (xl) an d the chiral susceptibility 
(Xm) \Karsch 02af . Both transitions would indeed be the same one or would take place at the same 
critical temperature. 
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The Stefan-Boltzmann law for bosons is [Landau 671 IGreiner 1)5] : 

eb = 3p = 9 Wf^~ () 

where g is the number of degrees of freedom due to spin, flavour, and colour charge of the 
considered particle. In consequence, e/T 4 or p/T A will be constant for such a matter. If 
one only considers photons (black body radiation) we obtain the Stefan-Boltzmann constant 
(g=2 for the two possible spins of the photon): 

= _^4 = 5m . 1Q -8 Wm -2 K -4_ (9) 

In natural units (temperature in MeV), the equation of a photon gas is 

e 7 = A x T 4 [MeV 4 ] (10) 

with A ~ 0.65. The Stefan-Boltzmann law for fermions is similar to that for bosons 
[Landau 67l IGreiner 95] : 

77T 2 {k B Tf 

e ' = 3p = ^^40- (1) 

The total energy density of this matter will be: e = e 7 + q + e g + e q ; with g = 8 for 
leptons (2 for spin, 2 for flavors and 2 particle-antiparticle), g = 16 for gluons (2 helicity 
states and 8 colour charges) and g = 36 for quarks (2 for spin, 3 colours, 3 flavours, and 2 
particle-antiparticle) : 

e={Ay + Ai + A g + A q ) x T 4 [MeV 4 ] (12) 

with A 7 =0.65, A h = 2.30, A 9 =5.26 and A q =10.36. 

For a small size plasma (radius below 10~ 10 m) or short lifetime, electromagnetic particles 
like photons and leptons could not reach thermalisation. They will be radiated by the 
thermalised medium but they will not be in equilibrium with the medium. Ignoring them, 
one gets e = 15.62 x T 4 [MeV 4 ] for 3 flavors of massless quarks and 8 gluons (see the value 
of e 5B /T 4 in Fig. §. 

F. Probes of the QGP 

1. Thermal radiation 

Thermal radiation from a QGP will allow to study several properties of the QGP like 
its temperature T. On the surface of the QGP volume, photonqj will escape. This is the 
thermal radiation. 

As we have estimated for an ideal QGP, the partial pressure of these photons on the QGP 
surface, will be given by the expression p = e/3 = 0.22T 4 and their energy distribution by 
the Planck law. In consequence the differential partial pressure dp/dE^ in natural units will 
be : 



dp(E 7 ,T)_ nnoA E, 



0.034 — ~' [MeV 3 ] (13) 

dE^ exp (£ 7 /T) - 1 [ J v ; 

where E^ is the photon energy in MeV. Considering massless particle and a QGP of a 
radius 7 fm and 10 fm/c lifetime, the thermal radiation spectrum is presented in Fig. M 



but also the other fundamental particles. 
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FIG. 9. Thermal photon production in a QGP of 7 frn radius and during 10 fm/c, as expected 
from a black-body radiation. 

for temperatures of 200 MeV, 500 MeV and 700 MeV. The corresponding photon yields for 
energies above 1 GeV are 52, 4600 and 16000, respectively. 

Obviously, the numerical example presented here is unrealistic since a 7 f m radius QGP 
will be transparent to photons. Under these conditions, the electromagnetic radiation of 
a thermalised QGP is not in thermal equilibrium with the medium which is producing it. 
Once a photon is produced, it will escape from the QGP, therefore the emission is from 
the volume and not from the surface as in the black-bo dy radiation. The c alculation of the 
thermal photon radiation from a QGP is complicated |Gelis 03| lArleo 03] . At first order, 
one could expect a reduction of the total number of photons emitted following the ratio of 
the strength of the strong and the electromagnetic forces o>qed/®-qcd- Only for large size 
QGP, with a radius above ~0.1 A, the black-body radiation model would become valid. 



2. Screening of the colour potential between heavy quarks in the QGP 



As we have already mentioned, the transition to the QGP only concerns the light quarks 
u, d and s, for which the chiral symmetry is a good approximation. Since heavy quarks 
explicitly break the chiral symmetry, they are not directly concerned by the transition to 
QGP. In other words, the bound states of heavy quarks (quarkonia) are not necessarily melt 
in a QGP and they could exist as bound states. For this reason, the se bound st ates become 
very interesting probes for measuring the temperature of the QGP |Matsui 86] . 

Let us see qualitatively which are the properties of a quarkonium embedded in a QGP. 

Quarkonia are bound states between two heavy quarks QQ: cc for the family r] c , J/i/j, 
ip(2S), Xc ■■■ and the states bb for the family T's and Xb- The bound state ti has not been 
experimentally observed and it will surely not exist due to the short lifetime of the top 
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quark. Finally, one should note that the hadronic states cb (cb) will have similar properties 
as quarkonium, although their decay should be similar to that of a B hadron. In vacuum, 
the quarkonium spectrum can be described via non relativistic models based on a potential 
interaction like: 

V{r) =ar-- (14) 

where a represents the string tension QQ and a is a Coulombian-like constant [Matsui 86J . 
For simplicity, let us assume that the potential is only Coulombian, so o = 0. 

If the QQ state is embedded in a QGP at a temperature T, the interaction potential 
between the heavy quarks will be affected by the presence of the free colour charges in the 
QGP. This is the screening of the potential. This phenomenon is well known in electromag- 
netic plasma. In the plasma, the Coulombian potential has to be replaced by a potential 
with a screening constant: 

V(r) = -- x e { - r/XD) (15) 

r 

where Xd is the Debye length. Let us assume that the average distance between the heavy 
quarks in a IS quarkonium state (J/ip or T(1S)) can be estimated by the Bohr radius 
expression: 

r B = — • (16) 

arriQ 

As a numeri cal example, one can consider for the J/ip m c =1250 GeV and a(m c ) = 0.36 
so r B = 0.44 fm. For the T{1S), m c =4200 GeV and a{m b ) = 0.22 



|Beringer 12 



|Beringer 12 



so tb = 0.22 fm. 



If tb -C Xd, the potential between the heavy quarks can be considered as a Coulombian 
potential and the bound state exhibits the same properties in the QGP as in the vacuum. 
However, if r# > Xd, the quarkonium properties will be modified by the medium, and it 
could happen that the quarkonium becomes an unstable state and therefore would melt. For 
electromagnetic plas mas, the De bye length depend on the temperature of the plasma and 
the charge density p [Stocker 99] : 

Ad = \^— (17) 

y 8nap 

Assuming that the previous expression is also valid for the QGPp| and an ideal ultra- 
relativistic gas p oc T 3 , one obtains : 

Ad ~ Tib- (18) 

And therefore, the quarkonium could be melt for temperature above T^: 

Td ~ f^WS ■ (19) 

A/87ra(T)rB 

For <x(T) ~ 0.2, one obtains that T d ~ 200 MeV (1.3T C ) for the J/-0 and T d ~ 400 MeV 
(2.6T C ) for the T(IS'). Assuming that for 2S states the r# would be twice larger, one would 
conclude that the dissociation temperature for *&' is ^ T c and for T(2S) similar to that of 

J/V. 

Of course, these calculations are qualitative and obtained without much detail, but they 
allow to show what is the role of parameters like the mass of the quarkonium and the 
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This assumption is not justified, but the conclusions that will be obtained are still valid. 



Bound state 
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Xc V J/V> T(25) ^ T(15) ' 



< T c < T c ~1.2T C ~1.2T C ~1.3T C ~2.0T C 



TABLE II. Upper bound of dissociation temperatures Tj of quarkonium states in units of the QGP 
transition temperature T c obtained by A. Mocsy and P. Petreczky in [Abreu 08] . 

strength of the interaction. This explains why ^/>(2S) resonance is easily melt with respect 
to J/ip and why T(1S) would melt at higher temperatures than that of J/ip. You will find a 
rigorous calculation of an upper bound of the dissociation te mperature o f quarkonium in the 
contribution of A. Mocsy and P. Petreczky in the reference |Abreu 08] . Their dissociation 
temperatures are quoted in table [TT} 

3. Parton - QGP interaction 

The QGP could also be studied via its tomography using high energy partons. QCD 
predicts that high energy partons will lose energy via gluon radiation when crossing the QGP. 
The order of magnitude of the parton energy-loss in QGP would be about AE ~1 GeV/fm 
and it is expected to be proportional to the gluon density. In addition QCD also predicts 
that the formation length of the radiated gluon will be larger than the average distance 
between the gluons in the QGP (interaction centres of the incident high energy parton). 
As a consequence several interaction centres will participate in the gluon emission from 
the parton, and the amplitude from the interaction centres will interfere (this phenomenon 
is called Landau-Migdal-Pomeranchuck effect) since the radiated gluon will be coherently 
emitted along all its formation length. For this reason for QGP thicknesses about 1-3 fm, 
the AE s hould be prop ortional to the square of the transversed path length in the QGP 
|Baier 97|IZakharov97] : 

AE~a s xC R x q( Pg ) x L 2 (20) 

where a s is the strength of the strong interaction, Cr is the colour charge factor q is the 
transport coefficient which depends on the gluon density (p g ) of the QGP and L is the 
thickness of the QGP. 

The energy lost will depend on the nature of the parton: 

• Gluons will exhibit larger energy-loss per unit of length than that of quarks. A relative 
factor 9/4 due to the colour charge, is associat ed to the g luonsstrahlung mechanism 
from a gluon with respect to that from a quark |Peigne 06] . 

• Heavy quarks are expected to lose less energy than light quarks, due to the absence 
of gluon radiatio n at forward an gles, below 9 < M/E, where M is the quark mass 
and E its energy [Dokshitzer Olj . This phenomenon, predicted by the QCD, is called 
dead-cone effect. The dead-cone effect should become measurable for beauty quarks, 
whereas this effect should remain relatively small for charm quarks. Moreover, elastic 
collisions with partons in the QGP could also contribute to the energy-loss of heavy 
quarks in the QGP. Finally the hadronization time scale for heavy quark hadroniza- 
tion increases due to its larger mass and it could occur, namely for the beauty, that 
hadronization takes place when the heavy quark is still traversing the QGP. 

One can wonder if other high-energy elementary particles like photons, electrons, elec- 
troweak bosons etc... could also be used to study the QGP. Photons and electrons will 
only interact electromagnetically and they should lose energy like in ordinary matter via 
bremsstrahlung emission and the production of electrons and positron pairs. However, the 
expected energy-loss is relatively small for QGP of a radius of tens of femtometers, about 
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~1% of the energy of the particle |Peigne 06|. In the case of electroweak bosons, they will 
decay quickly due to their short lifetime and only their daughter particles will interact with 
theQGP. 

III. HEAVY ION COLLISIONS AND HEAVY ION ACCELERATORS 

The study of hadronic matter in the laboratory is one of the challenges of experimental 
nuclear physics since the eighties. Today, the unique experimental method consists in accel- 
erating and colliding two heavy nuclei. In laboratories like CERN (Geneva, Switzerland), 
BNL (New York, USA), GSI (Darmstadt, Germany), and GANIL (Caen, France), nuclei 
are accelerated at energies that range from MeV to TeV beam energies. Depending on the 
center-of-mass energy of the collision, different domains of the phase diagram of hadronic 
matter can be studied. Before the collision, the nucleus-nucleus system is out of equilibrium. 
During the collision, the strong interaction between the constituents may dissipate a fraction 
of the available center of mass energy into the internal degrees of freedom of the system, and 
hopefully, a microscopic drop of hot hadronic matter could be created in the laboratory. The 
pressure gradient between the drop and the surrounding vacuum would be incredibly high 
and the drop will suffer a dramatic expansion against the vacuum. The temperature of the 
system will change during the expansion and a series of ephemeral thermodynamical states 
will be created. The complexity of this dynamical evolution of the system makes much 
more difficult the study of the intrinsic properties of the hadronic matter and a rigorous 
methodological approach has to be undertaken: 

• Collision dynamics. The systematic study of the different colliding systems, center 
of mass energies and impact parameter will be of vital importance; 

• Experimental probes. This implies that one can detect, identify and measure the 
kinematic properties of all the particles produced in the nucleus-nucleus collisions. 
This has not been always possible, and only large scale experiments in colliders like 
STAR, PHENIX, ATLAS, ALICE or CMS are able to perform such a complete mea- 
surement. This is the only way to measure all the experimental probes, like parti- 
cle multiplicity, light and strange hadron yields, transverse momentum and rapidity 
distributions, hadron correlations, azimuthal asymmetries, heavy quarks, quarkonia, 
direct photons, jets, dijets, electroweak bosons, photo-jet and electroweak bosons-jet 
correlations, etc ... 

• Experimental probes in cold nuclear matter. In addition, one has to study the 
experimental probes when the microscopic drop of hot hadronic matter is not created, 
namely, in peripheral collisions and/or induced-proton collisions; 

• Global interpretation. The results obtained have to be interpreted in one single 
scenario that explains coherently the whole phenomenology of experimental results. 

In order to create a drop of QGP in the laboratory, energy densities of about 1.0 GeV/fm 3 
have to be reached. Nucleus-nucleus collisions at relativistic energies have become a unique 
experimental method. Naively, one could assume that all the available energy in the center- 
of-mass is dissipated, during the collision, into the internal degrees of freedom of the nucleus- 
nucleus system. The latter statement is certainly a bad hypothesis since anyone will expect 
that a non negligible fraction of the available energy in the center-of-mass will not be dissi- 
pated to create hot matter. But this hypothesis allows to estimate the beam energy below 
which, the QGP cannot be formed. Under this hypothesis, the energy density of the drop 
will approximately be given by 

_ y/2E b x m x A 

V K } 
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where m is the mass of the nucleoli, Ey, is the beam energy in the reference system where 
one of the nucleus is at rest, A the atomic number and V the initial volume of the system. 
Assuming: 

V « 4/3tt x (1.124) 3 x A [fm 3 ], (22) 

we conclude that for beam energies Eb below ~20 GeV per nucleon (that is, a center-of-mass 
energy below a/snn=6 GeV per nucleon pair) the energy available in the center of mass is 
insufficient to heat a nucleus to energy densities above 1 GeV/fm 3 . It seems hardly possible 
that a drop of QGP could be then formed. Even at Eb ~ 20 GeV per nucleon, the stopping 
power of nuclear matter would not be strong enough to stop both nuclei, and in consequence 
only a fraction of the initial beam energy could be dissipated into the internal degrees of 
freedom of the system. Therefore Eb noticeably larger than 20 GeV per nucleon would be 
needed to reach the critical density of the QGP phase transition. 

At the beginning of the 80's, the American physicist J.D Bjorken imagined a scenario 
where the QGP would be efficiently formed. He described what would be the initial energy 
density and its evolution with time Bjorken 83] . As we will see later, one of the hypothesis 



of this scenario is only corroborated for Eb larger than ~250 GeV, that is an available energy 
in the center of mass larger that 25 GeV per nucleon pair. In addition, the QGP would be 
formed at baryonic potentials close to zero under this scenario. 

Therefore, it remains an open question whether the critical energy density could be 
reached or not, in the intermediate domain between £'6=20-250 GeV ( v /snn=6-25 GeV). As 
I will mention later, the results from the SPS experimental heavy ion program (1986-2000) 
at v /snn=17-19 GeV, hinted at the existence of a new state of matter in which quarks, 
instead of being bound up into more complex particles such as protons and neutrons, are 
liberated to roam freely. 

A. The Bjorken scenario of heavy ion collisions at ultra-relativistic energies 

At ultra-relativistic energies, nuclei are seen as pancakes in the center-of-mass system, 
due to the Lorentz contraction. The crossing time of the nuclei can be estimated as 

T cross = 2Rh, (23) 

where 7 is the Lorentz factor and R the radius of the nuclei. 
Bjorken assumed the following hypothesis: 

• The crossing time r cross is smaller than the time scale of the strong interaction. The 
latter can be estimated as T str0 ng ~ 1/Aqod ~ 1 fm/c. For a nucleus-nucleus collision, 
r cross is larger than r str0 ng only if 7 < 12. That is an energy in the center of mass 



above a/snn >25 GeV (so Ef > 250 GeV for a fixed-target experiment) 11 Under 
this hypothesis, the particles generated by the strong interaction between the 
partons, are created once the nuclei have already crossed each other. 



nucleon 



The distribution of the particle multiplicity as a function of the rapidity is assumed 
to be uniform. This is partially verified experimentally in Au-Au collisions at RHIC 
[PHOBO S 06|, in proton-antiproton collisions at Tevatron |CDF 90j . and at SPS ener- 
gies |Bjorken 83| . This condition ensures a rapidity symmetry of the system, allowing 
to create a uniform energy density in different rapidity slices, which simplifies consid- 
erably the description of the hydrodynamical evolution of the system. 



11 



To be noted that only RHIC and LHC colliders validate this hypothesis with 7 values of 100 and 1376 
(and 2750 after year 2015) respectively. 
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FIG. 10. Bjorken scenario \Bjorken 83\j for the formation of hot QCD matter. After a formation 
time Tf orm a volume with a high energy density is created. After equilibration at Tther, the evolution 
of the hot QCD matter follows the laws of the relativistic hydrodynamics. First, there is a longi- 
tudinal expansion until the system reaches a longitudinal size close to its transverse size, then a 
tridimensional expansion starts until the density is so low that no more inelastic (elastic) collision 
takes place. The system reaches then the so called chemical (kinetically) freeze-out. Finally all the 
particles will fly decaying to their daughter particles or reaching the detector. Typically only charged 
pions, charged kaons, protons, neutrons, photons, electrons and muons will reach the detectors. 
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Let's consider the volume centred in the nucleus crossing plane, at a time r after the 
nucleus crossing. This volume has a cylindrical shape with a thickness 2 Ad along the beam 
axis direction and a radius R ~ 1.124A 1 / 3 in the transverse plane. This volume will contain 
all the particles produced with a speed along the beam axis (z axis) below f} z < Ad/r. 
Since (3 Z = tanh (y) ~ y for y — > 0, the rapidity range Ay around y = of particles with a 
(3 Z < Ad/r will be 

Ay = . (24) 



and the total energy in the volume considered will be 



E 



dE 



dy 



2Ad 



x 



y=o 



T 



(25) 



where dE /dy is the total energy created by the strong interaction between the nuclei at y=0. 
For other rapidity domains, the previous expression can be easily generalised replacing the 
total ener gy b y the transverse energy E^. Finally, we can calculate the energy density in 
the volume^ 



<y) 



dE! 



dy 



TTiZV 



(26) 



which links the energy density with the transverse energy produced per unit of rapidity. 



1. Formation. 



The initial energy density can then be estimated assuming the time scale needed for the 
production of particles, as r form ~ r strong ~ 1 fm/c Pj 

Bjorken estimated the energy density for heavy ion collisions at beam energies of the 
SppS collider at CERN, that were ^/snn ~ 500 GeV per nucleon paiiPH and he obtained 
that the initial energy density were about 2-20 GeV/fm 3 , largely above the critical energy 
density to form the QGP. One can redo the exercise for heavy ion collisions at Tevatron 
energies (a/snn ~ 1.8 TeV |CDF 88tlCDF 90] ). and then the initial energy density would be 
4-30 GeV/fm 3 . 

Note that in Fig. 10 only hot matter created around mid-rapidity and its evolution is 
presented. Indeed one should keep in mind that the hot hadronic matter is created in the 
full rapidity range where the particle density is high enough to reach equilibrium. At RHIC 
energies, this is about 5 units of rapidity and at LHC energies about 8 units of rapidity. In 
the laboratory system, the hot matter slices at larger rapidities are indeed narrower due to 
the Lorentz contraction. 



2. Thermalisation. 



The particles produced inside the volume considered will interact. At these energy den- 
sities, and assuming a mean energy {E)=500 MeV, e/{E) ~ 8 — 60 particles per fm 3 will 
be reached. The average path length of particles inside the volume can be estimated as 
A ~ 0.02 — 0.12 fm, if one assumes an interaction cross-section of 10 mb. One could hope 



12 There is a factor 2 difference with respect to equation (3) in the original publication of Bjorken 
|Bjorken 83 . It is a known typo error in the original publication. 

13 Other estimates that provide smaller r s trong in the range 0.2-0.5 fm/c can be foreseen PHENIX 05a . 

14 The energy of the collider SppS is close to the available energies at RHIC: 200 GeV per nucleon pair 
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that the system will thermalise at a time r = Tther- Note that this is a strong assumption 
that must be i) validated by the experimental results and ii) supported by theoretical cal- 
culations. Experimental results seem to agree with the assumption of a fast thermalization 
of the system, but the theory has not been able to explain how thermal equilibrium could 
be reached in such a short time scale. This reminds a fundamental question to be answered 
and it is still a challenge for QCD theory to describe the first instants of the nucleus-nucleus 
collision at ultra-relativistic energies. In principle, the initial state of the nucleus-nucleus 
collision is characterised by the interaction of two high-density gluon clo uds. In th is respect, 
classical limits of the QCD theory (like the Colour Glass Condensate |Gelis 11] ) seem to 
be the best theoretical tool to study this problem. The typical Bjorken x of the two gluon 
clouds is (x) ~ 10~ 2 at RHIC and (a;) ~ 1(T 3 at the LHC 

3. Longitudinal expansion. 

At stages r > r ther the system should evolve like a fluid, following the laws of the rel- 
ativistic hydrodynamics. First a longitudinal expansion will take place since the pressure 
gradient in the beam direction will be larger than that in the transverse plane. It is expected 
that the energy densit y will evolve as e ~ l/r n with 1 < n < 4/3, which is obtained from 
the hydro dynamic law |Bjorken 83| 

^ = - e -±P. (27) 

dr r 

and for an ideal ultra-relativistic gas, this becomes e = 3p and thus n=4/3. The longitudinal 
expansion stays as a good approximation for stages r < Ti ong ~ R. 

4- 3D expansion and freeze- out phase. 

For stages r > Ti ong the system will evolve via a 3 dimensional expansion until the freeze- 
out stage is reached. At freeze-out, particle density is low enough to assume that particles 
do not interact, travel in the vacuum, can decay and finally reach the detector. Naively, the 
freeze-out will take place when the average path length of particles is similar to the size of 
the system A ~ R. For a cross-section of 10 mb, this corresponds to 0.15 particles per fm 3 
and therefore an energy density of 

e gd ~ 0.15 fm" 3 x 0.5 GeV ~ 0.075 GeV/fm 3 . (28) 

It is then expected that the freeze-out takes place as a hadron gas phase. Note that for 
a freeze-out temperature of T ge i = 150 MeV, one gets e/T 4 ~ 1.2, which fits pretty well 
with the prediction of lattice QCD calculations of Fig. [5] Finally it is worth mentioning 
that elastic cross-section is larger than inelastic one and one expects to observe two different 
freeze-out stages: chemical and kinetic freeze-out ones. 

B. Heavy ion accelerators and colliders 

Developments in heavy ions beams at ultra-relativistic energies have been performed in 
parallel as that at intermediate energies, since the main technical limitation was the ion 
source and the heavy-ion injection at low energies. The first heavy ion beams at relativistic 
energies where produced at AGS (BNL, USA) and at SPS (CERN, Switzerland) in the 80's. 
The energy in the centre of mass was 5 and 18 GeV per nucleon pair, respectively. The first 
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heavy ion collider was RHIC, built at BNL, which provided the first Au-Au collisions at 
a /snn=130 GeV in June 2000 and reached in 2001 its nominal energy of v / snn=200 GeV. 
Finally, LHC provided its first heavy ion collisions of Pb beam at a /snn=2760 GeV, a 14-fold 
increasing step with respect to RHIC, in November 2010 and hopefully this will turn into a 
28-fold factor (5500 GeV) from 2015 onwards. Today, RHIC and LHC are developing their 
heavy ion programs which are foreseen until 2025. 

1. The Alternating Gradient Synchrotron at BNL 

The AGS synchrotron was built in 1957 and allows the acceleration of high intensity 
proton beams at 33 GeV. Several Nobel prizes were obtained (1976, 1980 and 1989) linked 
to discoveries at AGS: J/i[) discovery in 1974, observation of the CP violation of the weak 
interaction in 1963 and the discovery of the muonic neutrino (1962). Since 1986, the AGS 
synchrotron has been used to accelerate Si ions at energies of 14 GeV per nucleon, after the 
construction of the beam line to inject heavy ions in AGS from the Tandem Van de Graaf 
(built in 1970). The Si beam from the tandem has an energy of 6.6 MeV per nucleon. The 
construction of the AGS booster in 1991 allowed to increase the AGS beam intensity and to 
accelerate heavier ions like Au up to 11 GeV per nucleon. Negative Au~ ions are extracted 
from the source and accelerated by the tandem to 1.17 MeV per nucleon and stripped to 
a beam of Au +32 . This beam is then injected in the AGS booster where the Au ions are 
accelerated to 90 MeV/nucleon. Finally the Au beam is stripped and injected into AGS 
where it is accelerated to the nominal energy of 11 GeV per nucleon. For 14 years, several 
fixed target heavy-ion experiments took place, like E866, E877, E891, E895, E896, E910, 
E917 to study the hadronic matter at high temperature^) 

2. The Super Proton Synchrotron at CERN 

The SPS was built in 1976, allowing for proton acceleration until 500 GeV. First, protons 
are accelerated by a linear accelerator called LINAC2, and then injected into the booster 
of the PS (Proton Synchrotron) and finally they are injected into the SPS to reach their 



nominal energy of 500 GeV 16 From 1986, the new electron-cyclotron resonance (ECR) ion 
source allowed the injection of multi charged heavy ions in the CERN accelerator system 
(LINAC3, PS booster, PS and SPS). The beam leaving from an ECR ion source containing 
a Pb plasma, has an energy of 2.5 KeV per nucleon with an ion charge Q = +27, and 
they are injected in the LINACS3 linear accelerator reaching a beam energy of 4.2 MeV 
per nucleon. Then the beam is stripped via a thin C layer 1 /zm thick, and becomes a 
+53 Pb beam, which is injected in the PS booster and PS accelerator, reaching an energy of 
4.25 GeV per nucleon. The Pb ions are then fully stripped in an aluminium layer of 1 mm 
thick, and they are injected in SPS to reach an energy of 158 GeV per nucleon. This beam 
is finally directed to the experimental fixed-target halls in SPS north area (NA) in France 
or SPS west area (WA) in Switzerland, where heavy ion collisions at a/snn take place. In 
addition to Pb, other ions have been also accelerated at SPS. At the beginning of the SPS 
heavy ion program, beams of O and S were accelerated at energies about 60 and 200 GeV 
per nucleon, and in the last days In ions were used for the NA60 experiment. During 20 
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Note that it is not clear the AGS could form deconfined matter since the initial energy density could be 
below 1 GeV/fm 3 . 

Initially the SPS was a proton accelerator. But SPS became a proton-antiproton collider with to the 
additional injection of antiproton beam. The latter was attainable thanks to the stochastic-cooling tech- 
nique in the SPS ring. The first collisions pp in SPS took place in 1981 at a center of mass energy of 520 
GeV. Two years later, the electroweak bosons were discovered by the UA1 and UA2 experiments. The 
stochastic-cooling and the discovery of the W, and Z bosons was awarded with the Nobel prize of physics 
in 1984. 
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years, many heavy ion experiments were built, installed and contributed to the SPS heavy- 
ion physics programme: WA80, WA93, WA98, WA85, WA94, WA97, NA57, Helios-2, NA44, 
CERES, Helios-3, NA35, NA49, NA36, NA52, NA38, NA50 et NA60. In 2000, the analysis 
and interpretation of the obtained experimental results was almost finished and a CERN 



press released was organisec 17 They announced that the physical results of the heavy ion 
fixed-target SPS experiment NA44, NA45, NA49, NA50, NA52, WA97 / NA57 and WA98 
hinted at the existence of a new state of matter in which quarks, instead of being bound up 
into more complex particles such as protons and neutrons, were liberated to roam freely. 



3. The Relativistic Heavy Ion Collider at BNL 

The first Au-Au collisions at 130 GeV per nucleon pair took place in June 2000 in RHIC at 
BNL (USA). It was the first collider ever built for heavy ions. AGS is the injector of RHIC, 
via a two Au beams at 9 GeV per nucleon which circulate in two different rings in opposite 
directions. In RHIC collider, 60 beam bunches in each ring are accelerated to the nominal 
energy of 100 GeV per nucleon and stored in two rings of 3.85 km perimeter length. The 
bunches of the two beams can collide in 4 interaction points along the RHIC ring, reaching 
nominal luminosities about 2 ■ 10 26 cm" 2 s _1 , that is a Au-Au collisions rate of 800 Hz. 
Recently, RHIC has been upgraded and is able to provide 5-10 times more instantaneous 
luminosity. In addition, RHIC collider allows to study collisions of polarised protons at 
500 GeV, and collisions of d-Au, Cu-Cu, Au-Au and U-U in the energy range 20-200 GeV 
per nucleon pair. Since 2000, the four experiments at RHIC: STAR, PHENIX, PHOBOS 
and BRAHMS have developed a high quality physics program, producing a huge amount of 
experimental results. Today only the two major experiments: PHENIX and STAR are still 
active and taking data. 



4- The Large Hadron Collider at CERN 

The LHC at CERN uses SPS as injector. SPS was upgraded to generate a Pb ions beam at 
177 GeV per nucleon, that are accelerated to a beam energy of 1.38 TeV. LHC provided the 
first Pb-Pb collisions at 2.76 TeV in November 2010, increasing by a factor 14 the centre-of- 
mass energy at RHIC. In November 2011 a new heavy-ion run took place at the same energy 
and the nominal instantaneous luminosity was reached, ~ 5 ■ 10 26 cm" 2 s -1 . It is expected 
that the nominal energy, 5.5 TeV per nucleon pair, will be reached after the long shutdown 
during 2013-2014. In principle the instantaneous luminosity at LHC and beam lifetime is 
limited by the huge cross-section of i) electromagnetic production of electron-positron pairs 
where the electron is captured by the Pb ions and ii) electromagnetic excitation of the Pb 
nucleus giant resonance, leading to neutron emission. Both processes are responsible for 
the Pb beam loss at LHC energy. LHC will be upgraded in 2018 to increase by a factor 
10 the instantaneous luminosity of the Pb-Pb collisions. At LHC, three of the four LHC 
experiments participate in the heavy ion program: ALICE, ATLAS and CMS. ALICE is the 
only LHC experiment devoted to the study of QGP. LHC will provide the first proton-Pb 
collisions at the beginning of 2013. 
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http://press.web.cern.ch/press/PressReleases/Releases2000/PR01.00EQuarkGluonMatter.html. 
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IV. SOME BASES ABOUT COLLISION CENTRALITY AND THE NUCLEAR 

MODIFICATION FACTOR 



In Fig. 10, the Bjorken scenario is presented for a central (zero impact parameter, b) 
collision is presented. Actually, collisions at any impact parameter between b = and 
b = R 1 + R 2 (the sum of the nuclear radius) could occur in the laboratory. It turns out 
that most of the collisions are indeed peripheral collisions, since the probability density 
is proportional to b. In experiments, the centrality of the collision can be estimated on 
an event-by-event basis via any observable C that monotonically varies with the impact 
parameter of the collision. The observable C can be the charged particle multiplicity or 
transverse energy in a given pseudo-rapidity interval, or energy at zero degree (at rapidities 
close to the beam rapidity), etc... Let us assume that i) f{C) represents the distribution 
of the observable C for a sample of non biased nucleus-nucleus collisions, that ii) C(b)>0 
and that hi) C(b=0) = 0. The centrality class n% of the most central collisions consists of 
nucleus-nucleus collisions where the observable C G (0, C n ) and 

f n C?l f(C)dC , N 

" = 100x JT7^ (29) 

The n% most central collisions are usually referred to as the centrality class 0-n%. Therefore 
the reaction class m%-n% (m < n) is defined by the collisions where the observable C G 

One of the experimental methods to quantify the nuclear medium effects in the production 
of a given observable (Ob) is the measurement of the nuclear modification factor (-Raa) i n 
nucleus-nucleus (A-A) collisions, defined as: 

yOb 



where (N co u) is the average number of binary nucleon-nucleon collisions 18 and Y®£ {Y^ 3 ) 
is the invariant yield of the observable Ob in A-A (pp) collisions at a given (same) center- 
of-mass energy. In the absence of nuclear matter effects, the nuclear modification factor 
should be equal to unity for experimental observables commonly called hard probes (large 
Pt particles, jets, heavy-flavour, etc). A similar factor R°&, measured in p-A collisions, is 
crucial in order to disentangle hot and cold nuclear matter effects in A-A collisions. 

V. BRIEF SUMMARY OF THE EXPERIMENTAL RESULTS AT RHIC AND 

AT THE LHC 

Due to a lack of time, I have not been able to complete satisfactorily these proceedings. 
For this reason, I am giving here a brief summary of the main results from RHIC (12 years 
of heavy ion programme) and from LHC (after the two first years of heavy ion programme). 

A. Initial energy density 

The multiplicity of charged particles d N ch /dr] was measured at RHIC [PHOBOS 00, 
IPHOBOS 021 IPHENIX 05aL IPHENIX 05b] (as well as the transverse energy [PHENIX Oil 

18 The average number of binary nucleon-nucleon collisions can be estimated by the product of the av- 
erage nuclear overlap function (of the nucleus- nucleus collision) and the inelastic proton-proton cross 
section [Miller 071 Id'Enterria 03] . 
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FIG. 11. Charged "particle pseudo-rapidity density per participant pair for central nucleus-nucleus 



collisions. 
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are superimposed on the diffractive pp collisions as 



a function of ^/snn heavy-ion and pp data, respectively. Figure 3 in reference I ALICE lOdlj . 



ST AR 04b] ). The most central Au-Au collisions at 200 GeV generate more than 600 charged 
particles per unit of pseudo-rapidity at mid-rapidity, which should correspond to about 900 
(charged and neutral) particles per unit of pseudo-rapidityFj Using the Bjorken model one 
can estimate that the initial energy density at mid-rapidity amounts to about 5-15 GeV/fm 3 . 
In addition the charged pa rticle multip licity remains constant within ~10% for 5 units of 
pseudo-rapidity (|r/| ^ 2.5) [PHOB OS llj . In the most central Pb-Pb collisions at 2.76 TeV, 
CMS has measu red that t ransverse energy at mid-rapidity is about ~2 TeV per unit of 
pseudo-rapidity [CMS 12b] . ALICE and ATLAS measured about ^1600 charged particles 



per unit of pseudo-rapidity [ALICE lOal IALICE Hal IATLAS 12b] (see Fig. |11). Consid- 
ering the increase of the mean hadron px at the LHC, the initial energy density at LHC 
is about three times larger than in Au-Au at RHIC top energy: 15-30 GeV/fm 3 . In both 
cases, the energy densities are several times larger than the critical energy density to form 
deconfined matter. Assuming that the system quickly equilibrates, the initial temperature 
could be estimated from the lattice QCD, assuming jx b ~ 0, (see Fig. [5]) as 



T 4 [MeV 4 



200 3 x 10 3 
12^5 



x e [GeV/fm 3 



(31) 



which gives an estimate of the initial temperature of 240-320 MeV at RHIC top energy and 
310-370 MeV at LHC 2.76 TeV energy. 



B. Equilibration 

Integrated hadron yields at mid-rapidities were studied in both energy domains. It is 
observed that the hadron yield ratios can be successfully described by a statistical model 

19 Grosso-modo one assumes that the pions are the most abundant produced particles. Only the tt° is 
neutral, so the total number of particles per unit of rapidity can be estimated as 600 x 3/2. 
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FIG. 12. Left: Comparison of thermal model predictions with RHIC data. Right: Thermal model 
fits to ALICE data on hadron production in central Pb-Pb collisions. From reference \ Andronic 10j . 



[Andronic 04j . In this model, the expanding hot system hadronizes statistically at the freeze- 
out, and therefore the hadron yields are given by the following expression: 



■ri; 



V 



9i 

2tt 2 



p 2 dp 



exp[(Ei ~ m)/T\ ± 1 



(32) 



with (+) for fermions and (-) for bosons, T is the temperature, iVj is the total number 
of hadrons of the species i, V the total volume of the system, g^ is the isospin and spin 
degeneration factor, Ei the total hadron energy and /z, the chemical potential. Considering 
zero total strangeness and isospin of the system, one can consider /ij = /j, b where //& is 
the baryonic chemical potential. Therefore only two parameters are needed to predict the 
hadron yield ratios: the freeze-out temperature and the baryonic potential. The analysis 
of hadron yield ratios allows to extract a similar freeze-out temperature of ~ 160 MeV at 
RHIC and at the LHC (see Fig. 12). The baryonic potential is //& ~20 MeV at RHIC 



and, as expected, a lower /if, at LHCTmdeed close to zero [Andronic 091 1 Andronic 12] . The 
value of the temperature at chemical freeze-out is indeed very close to th e phase transition 
temperature as predicted by lattice calculations presented in section IIP One should notice 
that, at LHC energies, proton and antiproton yields norma lised to the p ion yields exhibit 
an anomalous behaviour that has to be further investigated |Andronic 12] . 

The azimuthal distribution of particles in the plane perpendicular to the beam direction 
is an experimental observable which is also sensitive to the dynamics of the early stages of 
heavy-ion collisions. When nuclei collide at finite impact parameter (non-central collisions), 
the geometrical overlap region and therefore the initial matter distribution is anisotropic 
(almond shaped). If the matter is strongly interacting, this spati al asymmetry is converted 
via multiple collisions into an anisotropic momentum distribution |011itrault 93] . The second 
moment of the final state hadron azimuthal distribution with respect to the reaction plane 
is called the elliptic flow (f 2 ): 



E 



d 3 N 



d 2 N 



d 3 p 2n pTdprdy 



oo 

1 + ^{ 2 ^ C ° S W-^)]} 



n=l 



(33) 



where \1/r is the reaction plane, defined by the beam axis and the impact parameter. 
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FIG. 13. The elliptic flow as a function of the transverse momentum measured by PHENIX 
and STAR collaborations for hadrons, pions, kaons, protons and Lambda baryons. Figure 10 in 
reference ISTAR 05df . 



The elliptic flow has extensively be en studied at RHIC [STAR QlllSTAR 05a . PH0B0Sl)2l 
IPHENIX 09bllSTAR lOallSTAR 12a| (see Fig. pL and recently at LHC energies [ALICE lObt 
IATLAS 12al ICMS 12al IALICE llbllATLAS 12c| 7 Indeed the predictions from hydrodynam- 
ical models explain quite well most of the measurements of the elliptic flow of light hadrons 
at low pt (pt < 2 — 3 GeV). The elliptic flow measurements have been one of the major ob- 
servations at RHIC, evidencing that : i) the created matter equilibrates in an early stage of 
the collision, and then it evolves following the laws of the hydrodynamics; and ii) the formed 
matter behaves like a perfect fluid [PHENIX 05al ISTAR 05bl [PHOBOS 05l IBRAHMS 05] . 
Furthermore, several works (see for instance reference Nagle 10] ) managed to extract values 
of transport properties, like the ratio of the shear viscosity over entropy from the exper- 
imental results. The conclusion was that the hot matter behaves as a perfect fluid and 
the mean free path of the constituents is clos e to the quan tum limit. ALICE presented 
the first elliptic flow measurement at the LHC [ALICE 10b] in agreement with other LHC 
results [ATLAS 12a] . It was observed a similarity between RHIC and the LHC of px- 
differential elliptic flow at low px, which is consistent with predictions of hydrodynamic 
models (px < 2 — 3 GeV/c). The elliptic flow is now being studied in much more details 
at LHC, for identified hadrons and as a function of the pseudo-rapidity (see for instance 
[AT LAS 12cj ). Preliminary results on these topics are intriguing and it is certainly too 
early to conclude about their interpretation. In addition, the ellipt ic flow is b eing studied 
for the first time at very high px, (up to px ~ 40 — 60 GeV/c) [CMS 12a] . Other mo- 
ments of the hadron azim uthal distribution have been studied at RHIC and at the LHC 
[ALICE lib} IATLAS 12aj . In particular non-zero v$, was observed and it results from the 
fluctuations of the initial spatial distribution of the energy density. Complementarily to 
the V2, the measurement of higher harmonics of the azimuthal distribution is also crucial to 
constrain the shear viscosity over entropy ratio in models. 
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FIG. 14. Invariant cross section (pp) and invariant yield (Au-Au) of direct photons as a function of 
pp. The three curves on the pp data represent NLO pQCD calculations, and the dashed curves show 
a modified power-law fit to the pp data, scaled by Taa- The dashed (black) curves are exponential 
plus the Taa scaled pp fit. The dotted (red) curve near the 0-20% centrality data is a theory 
calculation. Figure 3 in reference [PHENIX 10}. 



C. Initial temperature 



As we have seen in section II F 1[ if QGP drop is formed, it should emit thermal radia- 
tion in the high energy 7 domain. PHENIX collaboration have measured e + e~ pairs with 
invariant masses below 300 MeV/c 2 and 1<Pt < 5 GeV/c in Au-Au collisions at 200 GeV 
[PHENI X 10| . The most central Au-Au collisions show a large excess of the dielectron yield 



see Fig. 14). By treating the excess as internal conversion of direct photons, the direct 



photon yield is deduced. The yield cannot be explained by Glauber scaled NLO pQCD 
calculations. However, hydrodynamical models with an initial temperature of 300-600 MeV 
are in qualitative agreement with the data. The evidence for the production of thermal di- 
rect photons, with an initial temperature source above the QGP transition temperature 
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FIG. 15. Dimuon invariant-mass distributions in Pb-Pb (left) and pp (right) data at ^/snn= 2.76 
TeV. The solid (signal + background) and dashed (background- only) curves show the results of the 
simultaneous fit to the two datasets. Figure 1 in reference J CMS 123/ . 



represent an important experimental observation. Preliminary results from AL ICE abou t 
thermal photon production in central Pb-Pb at 2.76 TeV are already available [Wilde 12] . 
The supposed thermal photon yield exhibits a 40% larger inverse slope at LHC than that 
at RHIC. The latter is in qualitatively good agreement with the expected relative increase 
of the initial temperature from RHIC to LHC energies. 

Quarkonium was proposed as a probe of the QCD matter formed in relativistic heavy-ion 
collisions more than two decades ago. A familiar prediction, quarkonium suppression due 
to colour-screening of the heavy-quark potential in deconfmed QCD matter [Matsui 86] . has 
been experimentally searched for at the SPS and RHIC heavy-ion facilities. 

CMS collaboration has performed the first measurement of the upsilon resonances (T(1S), 
T(2S) and T(3S)) at the LHC [CMS llal |CMS 12g[ ICMST2c] . The results indicate a sig- 
nificant decrease of the T(2S) and T(3S) Raa (see Fig. 15). The T(1S) Raa is about 0.41 
for the most central collisions. One should note that about~50% of the upsilon production in 
hadronic collisions is ex pected to result from the radiative decays of higher bottomonium res- 
onances |Bedjidian 04| . If one assumes that high resonances are dissolved, one would expect 
to measure a nuclear modification factor for the T(1S) about 0.5. The present measure- 
ment would be compatible with a formation of a QGP at the LHC at an initial temperature 
between 1.2-2.0 times the critical temperatures (see Tab. [H]), so absolute temperatures be- 
tween 200-400 MeV. Since the melting temperature of T(2S) and 3 /if are expected to be 
similar (see Tab. [H]) one should expect a similar decrease of the 3 /if) Raa at LHC energies. 

The PHENIX experiment at RHIC reported the observation of J /if suppression in central 
Au-Au collisions at v /snn=200 GeV (10 times higher than the maximum energy in the CM 
at SPS) [PHENIX 07l IPHENIX llbl IPHENIX 12b] . Deuteron-gold c ollisions have been 
used to constrain cold nuclear matter (CNM) effects at RHIC energies [PHEN IX llaj . As 
a consequence, J/if) suppression due to dissociation in QGP matter is roughly estimated to 
be 40-80% in central Au-Au collisions at RHIC energies. Since about 40% of the J /if) yield 
results from the decays of higher resonances, it remains an open question whether the J /if) 
is melt or not melt at RHIC energies. Finally, the STAR experiment has measur ed a smaller 
suppression at high transverse momentum (p T > 5 GeV/c) at mid-rapidity [STAR 09al 
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FIG. 16. Left: ALICE px dependence of the inclusive J/ip Raa measured in Pb-Pb collisions 
at ^/snn = 2.76 TeV is compared to PHENIX and CMS measurements. Figure 2 in reference 
[Sui re 12f . Right: J/ip elliptic flow as a function of px in the centrality range 20-60% and in the 
rapidity range 2.5 < y < 4 (red squares). ALICE data are compared to STAR measurement per- 
formed in Au-Au collisions at ^/snn=200 GeV and in the rapidity range \y\ < 1 (black diamonds). 
ALICE data are compared with a transport model predictions. Figure 2 in IMassacrier 12f . 

ISTAR 12c] although the experimental errors remain large. 



D. The phase of deconfinement 

At LHC energies, on average one J /if} particle is expected to be produced in every cen- 
tral Pb-Pb c ollision, together with about 50-100 cc quark pairs. As suggested in 1988 
|Svetitsky 88 , under these conditions the charm quark yield per unit of rapidity could be 
large enough to enhance the charmonium production in later phases of the hot QCD-matter 
dynamical evolution, in particular when the energy density is low enough to enable the char- 



monium bound state to be formed IBraun-Munzinger 00[ Grandchamp 04 lAndronic 11] , 
The ALICE collaborat ion reported the first measurement of the J/ip nuclear modification 
factor at LHC energies [ALICE 12dj . Contrary to the expectations from T(25) suppression, 
J/i/j Raa was found to be about 0.5 in the most central Pb-Pb collisions and does not exhibit 
a significant centrality dependence. In addition, the J /if} Raa was found to be larger than 
that measured at RHIC. Contrary t o RHIC o bservation, the J/ip Raa is large at low px and 
then decreases with increasing px [Suire 12] (see Fig. 16 left). Finally, a hint of non-zero 
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elliptic flow was also measured by the ALICE collaboration [Massacrier 12] (see Fig. 
right). These experimental observations suggest that J/ip production at LHC energies is 
governed, for an important part, by charm quark recombination processes. The production 
of charmonium in latter stages of the QGP evolution would certainly be a direct probe of 
the deconfinement phase. 

Heavy-flavour hadrons, containing charm and beauty heavy quarks, are effective probes of 
the QGP. During the deconfined phase, low px heavy quarks will inter act with the medium 
modifying its initial kinematical properties [van Hees 08 [ IGossiaux 09] and, in the extreme 
scenario, they could become fully thermalised. The first ALICE results on the nuclear 
modification factor Raa for charm hadrons in Pb-Pb collisions indicate strong in-medium 
energy loss for charm quarks. The D°, D + , and D* + Raa, were measured for the first time 
as a function of transverse momentum and centr ality. The suppressio n is almost as large as 
that observed for charged particles |ALICE 12g[ [Conesa del Valle 12] (see Fig. 17). A hint 
of non-zero flow of D hadrons was also measured [Ortona i~2] . High precision measurements 
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FIG. 17. Average Raa of D mesons in the 0-20% centrality class compared to the nuclear 

modification factors of charged particles and non-prompt J/ip from B decays measured by the CMS 
collaboration in the same centrality class. Figure 8 right in reference 'ALIC E 12g h 



of heavy-flavour hadrons at low pt will remain an experimental challenge during the next 
10 years at RHIC and the LHC. 



E. The opacity of the QGP 

Heavy ion collisions at LHC energies allowed to study for the first time the interaction 
between hard partons produced in the first stage of the hadronic collisions, with the QGP. 
The easiest experimental way to address this topic was via the study of high px particle 
yields and high px hadron-hadron correlations. High pt particles are produced by the 
fragmentation of partons (quarks or gluons) in a time scale around r frag ps _E/Aq CD x 
R, where E is the energy of the parton and R is the typical size of a hadron (~ 1 fm). 
For energies above 5 GeV, the fragmentation time scale is about 20 fm/c. Therefore, in 
heavy ion collisions, partons are expected to fragment after traversing the QGP. One of the 
major discoveries at RHIC was the suppression of high p T hadron Raa and the quenching 
of back-to-back hadron co rrelations [STAR 021 IPHENIX 03a] IPHENIX 03b] IBRAHMS 04l 
STAR 04a, PHENI X 05a] . This observation has been explained by the formation of a QGP 
drop where the init ial har d partons interact losing a non negligible fraction of their initial 
energy (see section II F 3). QCD inspired mo dels assuming gluon radia tive energy loss in 



the QGP are in good agreement with the data [STAR 05b] IPHENIX 05a] . On this topic the 
phenomenology is very rich and m any experimental detailed studies have been performed 
[STAR lib! TSTAR 12dl ISTAR 12bl [PHENIX 081 IPHENIX 12aj . 

The first LHC results on Raa have confirmed RHIC results and extended the px ranges 
until values as high as 100 GeV/c (ALICE lid ICMS 12ellALICE 12a] . The results indicate 
a strong suppression of charged particle production in Pb-Pb collisions and a characteristic 
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FIG. 18. Left: Nuclear modification factor R\a(pt) for it in central (closed circles) and peripheral 
(open circles) Au-Au at ^/s NN = 200 GeV. Figure 3 from reference JPHENIX 03b)/ . Right: Raa of 
charged particles measured by ALICE in the most central Pb-Pb collisions (0-5%) in comparison 
to results from CMS and model calculations. Figure 4 from reference J ALICE 12a\j 



centrality and px dependence of the nuclear modification factors. In the most central col- 
lisions, the Raa is strongly suppressed (Raa ~ 0.13) at px = 6-7 GeV/c. Above p T = 7 
GeV/c, there is a significant rise in the nuclear modification factor, which reaches Raa ~ 
0.4 for px > 30 GeV/c (see Fig. 18). The latter is in good agreement with models based on 
radiative energy loss of gluons hiQGP. 

At LHC the studies of jets in heavy ion collisions becomes possible. The ATLAS col- 
laboration presented the first results on jet reconstruction in Pb-Pb collisions at the LHC 
[ATLAS TO] . Jets were reconstructed up to transverse energies of 100 GeV. An asymmetry, 
increasing with centrality, was observed between the transverse energies of the leading and 
second jets (see Fig. 19). This is an outstanding confirmation of the strong jet energy loss 
in a hot, dense medium, as it was inferred from the studies of the high p T Raa and hadron 
correlations at RHIC. Simi lar conclus ions were obtained from the measurement performed 
by the CMS collaboration [CMS lib] . 

At LHC, the phenomenology on studies related to QCD energy loss is also very rich. 
Many measu rements that are not des cribed here have been perfor med, like h adron-hadron 
correlations [ALICE 12ej . single jets [C MS 12hj and gamma-jets [CMS 12dj . In the next 
10 years, high precision measurements will be performed on these channels and other more 
exotic ones, like Z-jet, will be studied. 

The study of high px Raa of heavy flavour hadrons should shed light on the QCD energy 
loss mechanisms. According to QCD, the radiative energy lo ss of gluons sho uld be larger 
than that of quarks. In addition, due to the dead cone effect [Dokshitzer Olj . heavy quark 
energy loss should be further reduced with respect to that of light quarks. Many studies 
were performed at RHIC, mainly via the semileptonic decay of heavy flavour hadrons. A 
strong suppression was observed but quantitative conclusions are not yet avail able. At the 
LHC, ALICE collaboration has measured the high p T R A a of D°, D + , and D*+ JAL1CE 12g 



IConesa del Valle 12| and the high px Raa of semi-muonic decay of heavy-flavours (charm and 
beauty) [ALICE 12fj . The CMS collaboration has measured the high pt Raa of J /if) from 
beauty hadron decays. These results indicate strong in-medium energy loss for charm and 
beauty quarks, increasing towards the most central collisions. It seems that J /if) from beauty 
hadron decays are less suppressed than charm hadrons, but systematic uncertainties are still 
large. In the next 10 years, thanks to the upgrades of the LHC and RHIC experiments, 
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FIG. 19. Top: dijet asymmetry distributions for data (points) and unquenched HIJING with 
superimposed PYTHIA dijets (solid yellow histograms), as a function of collision centrality (left 
to right from peripheral to central events). Proton-proton data from -y/s = 7 TeV, analyzed with 
the same jet selection, is shown as open circles. Bottom: distribution of A<p, the azimuthal angle 
between the two jets, for data and HIJING+PYTHIA, also as a function of centrality. Figure 3 
from reference ] ATLAS 10j . 

higher precision measurements will become available. 



F. Other interesting measurements 

Among the huge amount of experimental results that have not been described in this 
section, I would like to quickly mention the following ones: 



The measurement of electro- weak boson Raa, proposed by |Conesa del Valle"08] . has 
become possible at LHC . CMS and ATLAS collaboration h as performed the first mea- 
surements at the LHC [CMS llcl IATLAS 12dl ICMS 12f| . These have been funda- 
mental measurements and (unfortunately) the measured nuclear modification factor is 
compatible with unity, as it was expected. 

The charged particle multiplicities measured in high-multiplicity pp collisions at LHC 
energies reach values that are of the same order as those measured in heavy-ion col- 
lisions at lower energies (e.g. they are well above the ones observed at RHIC for 
peripheral Cu-Cu collisions at 200 GeV [PHOBOS llj). Therefore, it is a valid ques- 
tion whether pp collisions also exhibit any kind of collective behaviour as seen in 
these heavy-ion collisions. An indication for this might be the observation of long 
range, near-side angular correlations (ridge) in pp collisions at 0.9, 2.36 and 7 TeV 
with charged particle multiplicities above four times the mean multiplicity |CMS 10] . 
Recently J/ip yields were measured for the first time in pp collisions as a function of 
the charged particle multiplicity density [ALICE 12c] . The study of high multiplicity 
pp and p-A collisions will be an exciting topic in the next years. 

Antimatter can efficiently be created in heavy ion collisions. STAR collaboration 



reported the first observation of the anti-helium-4 nucleus [STAR Hal 
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• Finally, ultra-peripheral heavy ion collisions at RHIC and at the LHC have be- 
come a powerful high lumi nosity phot on beam. Many interesti ng measurem ents of 
vector meson p |STAR 081 ISTAR 09bl ISTAR 11a] . multi-pions [STAR 10b] or J/<0 
jPHENIX 09al IALICE 12b] are being performed in both colliders. 

G. Caveat on cold nuclear matter effects 

This topic has not been addressed in the present proceedings. The study of cold nuclear 
matter effects in proton or deuteron induced collisions is of outstanding importance. Many 
of the interpretations of the experimental results given above can only be confirmed via 
the study of these collisions. At RHIC energies, deuteron induced collisions have been 
extensively studied. At LHC, the first run p-Pb has taken place beginning of 2013. 

VI. OTHER LECTURES ON QGP 

The following references that will certainly complement the present lectures: 

• Lectures of Larry MacLerran, The Quark Gluon Plasma and The Color Glass Con- 
densate: 4 Lectures [Mc.Lerran 01] . 

• Lectures of Frithjof Karsch, Lattice Results on QCD Thermodynamics [Karsch 02a] . 

• Lectures of Jean-Paul Blaizot, Theory of the Quark-Gluon Plasma [Blaizot 02] . 

• Lectures of Ulrich W. Heinz, Concepts of Heavy-Ion Physics [Heinz 04] . 

• Lectures of Anton Andronic and Peter Braun -Munzinger, U ltra relativistic nucleus- 
nucleus collisions and the quark- gluon plasma, |Andronic 04] . 

• Lectures of Thomas Schaefer, Phase of QCD [Schaefer 05] . 

• Lectures of Bernt Muller, From Quark-Gluon Plasma to the Perfect Liquid |Miiller 0"7] . 

• Lectures of Je an- Yves Ollitrault, Relativistic hydrodynamics for heavy-ion collisions 
jUllitrault 07] . 

• Lectures of Tetsufumi Hirano, Naomi van der Kolk and Ante Bilandzic, Hydrodynamics 
and Flow [Hirano 08] . 

Lectures of Carlos Salgado. Lectures on high- energy heavy-ion collisions at the LHC 
Salgado 09| . 



Article by Michael L. Miller, Klaus Reygers, Steph en J. Sand ers, Peter Steinberg, 
Glauber Modeling in High Energy Nuclear Collisions [Miller 07] . 

Article by David d'Enterria, Hard scatteri ng cross sectio ns at LHC in the Glauber 
approach: from pp to pA and A A collisions [d'Enterria 03] . 

Article by Berndt Muller, Jur gen Schukr aft and Bolek Wyslouch, First Results from 
Pb+Pb collisions at the LHC [Muller 12] . 



In French: Proceedings of Joliot-Curie School in 1998: http : //www . cenbg . in2p3 . f r/ 
heberge/EcoleJoliotCurie/coursJC/J0LI0T-CURIE%201998.pdf [Joliot-Curie 98] . 
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In French: Proceedings of Joliot-Curie School in 2005: http : //www . cenbg . in2p3 . f r/ 
heberge/Ec ole JoliotCurie/cours JC/ J0LI0T-CU RIE%202005 . pdf [Joliot- Curie 05 



• In French: My HDR (Habilitation a Diriger des Recherches) [Martinez 06j . 
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